Universality of the critical conductance distribution in various dimensions 
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We study numerically the metal - insulator transition in 
the Anderson model on various lattices with dimension 
2 < d < 4 (bifractals and Euclidian lattices). The criti- 
cal exponent v and the critical conductance distribution 
are calculated. We confirm that v depends only on the 
spectral dimension. The other parameters - critical dis- 
order, critical conductance distribution and conductance 
cummulants - depend also on lattice topology. Thus only 
qualitative comparison with theoretical formulae for di- 
mension dependence of the cummulants is possible. 

PACS numbers: 71.30.+h, 71.23.An, 72.15.Rn 

It is commonly accepted, though not proved, that metal 
- insulator transitions (MIT) can be described by one- 
parameter scaling theory jjj. The critical exponent v 
which describes the divergence of correlation length at 
MIT depends only on the system dimension for a cho- 
sen universality class. Microscopic details of models do 
not affect it. This was confirmed by numerical analysis of 
quasi-one-dimensional (Q1D) systems M. Theoretical de- 
pendence of v on dimension d = 2 + e was found in |, | 
for small e. Numerically, v(e) was studied on bifractals 
I- 

The conductance g was originally chosen as the order 
parameter in the scaling theory Jl]]. Soon it became clear, 
that the absence of self-averaging of g in the critical re- 
gion must be taken into account || f?J. The shape of the 
critical conductance distribution P(g) in 3D models was 
numerically analysed in detail || ^|, |l0|, |ll|, [l^]. Contrary 
to the critical exponent, P(g) is not universal. Its shape 
depends not merely on the dimension JlO[ ] and physical 
symmetry but also on boundary conditions [ |i"2| and 
even anisotropy p3| . Nevertheless, for a given physical 
model the mean conductance and resistance follow one 
parameter scaling 

Analytical theory of MIT is restricted to systems with 
dimension close to the lower critical dimension: 2 + e with 
£ <C 1 jl5| . In spite of predicted non-universality of higher 
order conductance cummulants {Sg n ), 



the distribution Pig) should be universal in the infinite 
system size limit jl6|. For small e, the bulk of the dis- 
tribution is approximately Gaussian near the mean value 
(g). The parameters of Gaussian peak, 



(9) 



and 



varg 



{9? 



(2) 



agree with the estimation of the first cummulants (Q). 
For large g, theory predicts long power-law tail P(g) ~ 
<7 _1_2 / £ |L6|. Numerical results for 3D systems show com- 
pletely different P(g), indicating that theory is not appli- 
cable to large e J§|]. 

In this Letter we present the critical exponent and the 
critical conductance distribution for Anderson model ob- 
tained numerically on three bifractal lattices with dimen- 
sion 2 < d < 3. They all possess the same fractal dimen- 
sion df = log 3/ log 2 + 1. Two of these lattices have the 
same spectral dimension d s . Their critical exponents are 
identical within error bars. This confirms the universality 
of MIT. The shape of the P(g), and the value of the criti- 
cal disorder, depend not only on d s but also on the lattice 
topology. This novel non-universality of P(g) is found also 
by numerical analysis of two different 3D lattices. 

Topology dependence of the P(g) and of the critical 
disorder disables the verification of theoretical formulae 
for conductance cummulants. Known theoretical formulae 
may be valid only for d-dimensional hyper-cubes which can 
be numerically simulated only for integer d. Thus we cal- 
culated the critical parameters for four dimensional (4D) 
lattice and compared them with known results for 3D lat- 
tice. Ratios of the first two conductance cummulants are 
(<?)3d/ (g)4D ~ 2 and var g^ ex var (74D, in full agreement 
with relations (0). 



We consider the Anderson Hamiltonian 

H = J2s n \n)(n\+J2\n)(n'\ + \n'){n\ 

n [tin'] 



(3) 



(Sg n ) = 



-n-2 



n < no — e 
n > e^ 1 



(1) 



Random energies e n are uniformly distributed from 
(— W/2, W/2). Parameter W measures the disorder 
strength. Fermi energy equals to zero, n numbers lattice 
sites, and \nn'] are two nearest-neighbor sites. 

All systems under consideration are linear in the z- 
direction. In the plane perpendicular to the current di- 
rection we construct fractals A-C according to figure 1 . A 
combination of a c?-dimensional fractal with linear chain 
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in z direction produces a bifractal with dimension d + 1 
both for d = d s and d — df. The length of the system in z 
direction is L = 2™ and the number of lattice sites in the 
slice grows as N — 3™ for the nth generation of bifractal. 
We study also two 3D models with triangular and honey- 
comb 2D lattice in the xy plane (referred as 3t and 3h) 
and compare the obtained results with known data for 3D 
cubes (3s) 0, g |14 



In U it was supposed that critical parameters are com- 
pletely determined by the spectral dimension d s of the 
lattice Jl|. Following (l| and 0] we find the analyti- 
cal values for the spectral dimension for fractals A-C po| : 
df = df = 2 log 3/ log 5 and df = 2 log 3/ log 6. The last 
value differs slightly from the one obtained by numerical 
simulation of a random walker [^J. 

Linear form of all systems in z direction enables us to 
apply the standard numerical procedure for calculation of 
the critical disorder W c and critical exponents v in Q1D 
systems and the Landauer formula for conductance (in 
units e 2 /h) H 

g = 2Tr tH. (4) 
The same spectral dimension of fractals A and B re- 



quires v ■ 
2.2 ± 0.2 (17 



v . Although recently reported data, 



and 



2.5 ±0.25 1 



could be regarded 
as equal to each other within error bars, we wanted to 
check their equivalence more accurately. In our analysis 
we consider first five generations (L < 32) of fractals. The 
smallest Lyapunov exponent Z\ was calculated with accu- 
racy of 0.1% for L=4,8 and 16 and 0.5% for L = 32. Our 



result indeed confirms v J 



(Table 1). 



We tested the universality of v also for 3D lattices with 
various topology. Q1D systems up to 14 2 lattice sites in 
the xy plane were considered. The accuracy of the first 
Lyapunov exponent was 0.1% for small crossections, and 
decreases as L increases, being only 0.5-1% for L = 14. 
Our results confirm the universality of v as expected. 

For 4D hyper-cubes we involve systems up to 7 3 x oo. 
Our resulting v rs 1.1 is in a very good agreement with 
the previous ones ^ ^2j . The obtained critical parameters 
are presented in Table 2. 

The critical conductance distribution P(g) for lattices 
A, B and C is presented in figure 2. For bifractals A and 
B we can approve the system-size independence of P(g). 
For C, we do not reach the limiting form of P(g) because 
of finite-size effects which are much stronger at lower d s . 

As supposed, the mean conductance (g) increases as d s 
decreases. We cannot, however, compare our data with the 
theory, since the theoretical analysis has been performed 
only for (i-dimensional hyper-cubes. We can only describe 
some general features of the critical distribution: (i) P{g) 
converges to Gaussian as d s decreases to 2, as predicted in 
PI . (ii) We found no evidence of the power law decrease 
of P(g) for g 3> 1. This could be caused by the small 
statistical ensembles (we have only ~ 1.000 samples in the 
6th generation, L = 64). We note that the qualitative 
arguments against the power-law decrease [0 may not be 



valid in the limit d — > 2 + because the differences between 
Lyapunov exponents Zi+i — Zi become very small for e — > 0. 

Different forms of P(g) for bifractals A and B indicate 
that the critical conductance distribution and also con- 
ductance cummulants depend on the lattice topology. A 
more convincing proof of this statement is in figure 3 where 
we present the critical conductance distributions for 3D 
systems with honeycomb and triangular lattice in the xy 
plane and compare them with P(g) for 3D cubic lattice. 

Figure 4 presents P(logg) for 4D cubic lattice (with 
fixed boundary conditions) . The mean conductance is half 
of that in 3D cubic lattice. The shape of P(g) and P(log g) 
is similar to the one for 3D (figure 3) and it can be analyzed 
by standard methods jl0[ ^3) . 

Finally, figure 5 summarizes the dimension dependence 
of the critical parameters. The main difference between 
our numerical data for v and those published in || is the 
estimation of the spectral dimension which is more accu- 
rate in the present work. Contrary to Ref. |^], we do 
not try to fit the dimension-dependence of our data to 
any simple function. As pointed out in e-dependence 
of v is non-trivial. The leading term v — 1/e would be 
observable only on lattices with much smaller spectral di- 
mension. As an exact knowledge of d s is crucial for nu- 
merical estimation of the e-dependence of v 1 we did not 
analyse statistical fractals, for which d s must be estimated 
by numerical simulations. Another disadvantage of statis- 
tical fractals is that due to the topology dependence of 
P{g) we do not suppose that system-size invariant critical 
conductance distribution can be found on available scales. 

The dependence of the conductance cummulants on the 
lattice topology disables quantitative comparison of nu- 
merical and theoretical data. We can only conclude that 
the mean conductance increases as e decreases. To our 
surprise, (g) and var g for 3D and 4D seem to follow (Q). 
Does it mean that the relation no = 1/e in ([!]) under- 
estimates the upper bound no? A positive answer could 
explain the absence of the power law decrease of P(g) for 
large g. 

In conclusion, we have analyzed the critical parameters 
of the metal-insulator transition on lattices with different 
dimension and different topology. We confirm that the 
critical exponent v depends only on the spectral dimen- 
sion of lattice. This confirms universality of the MIT. Our 
result for v in 4D systems agrees with the previous nu- 
merical estimations H [22], but differs considerably from 
theoretical expectation v = 1/2 pij . 

We present the critical conductance distribution on the 
lattices of dimension 2.226 < d s < 4 and compare them 
with theory. We prove that two lattices with the same 
spectral dimension but different lattice topology have a 
different critical conductance distribution. This prevents 
a quantitative comparison of numerical data for the con- 
ductance cummulants and for the critical disorder with 
theoretical formulae. In agreement with Jl6|], P(g) con- 
verges to the Gaussian when dimension decreases toward 



2 



the lower critical dimension d s = 2. We found no power- 
law tail of the distribution, maybe due to the restricted 
size of the statistical ensembles. 

We thank L. Samaj for fruitful discussions. This work 
was supported by Slovak Grant Agency VEGA, Grant n. 
2/7174/20 
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Figure 2: Critical conductance distribution P(g) for 
bifractals A, B and C. The number of samples in statis- 
tical ensembles decreases from 10 5 for smaller systems to 
~ 1.000 for L — 64. The last figure shows the system size 
dependence of mean conductance (open symbols) and var 
g (full symbols). Note that finite-size effects increase as 
critical disorder decreases. 
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Figure 4: Critical conductance distribution of log g for 
4D hyper-cubes with fixed boundary conditions. The pa- 
rameters of the statistical ensembles are presented in the 
legend: L, (logg), var log g and number of samples in 
the statistical ensemble. Inset: The first Lyapunov ex- 
ponent z\ as a function of disorder for the Q1D lattice 
L 3 x L z and 4 < L < 7. Solid lines are linear fits 
z 1 (W,L) = z[ 0) (L) + Wz[ 1 \l) 



10" 



D)icr 
o 



10" 



10 





Triangular 




,6 






r / 


a. 


\ < 






2 










a 




□ 






<> 






* / i 


Honeycomb 








o° a 
, s 1 * 




t 






□ O 




i 




° / ^ 


A L=10 






A 








O 










□ 


4 n 









-10 



-6 -4 

logg 



Figure 3: Critical conductance distribution for three di- 
mensional systems with triangular and honeycomb struc- 
ture in the xy plane. For comparison, P(g) for the cubic 
lattice is also presented (dotted line). 
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Figure 5: Dimension dependence (e = d — 2) of critical 
exponent (left) and of the mean and variance of conduc- 
tance (right). Systems with the same spectral dimension 
d s have different mean and variance of conductance, but 
they possess the same critical exponent v. For details see 
Tables 1 and 2. 
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bifractal ABC 



d s 


2.365 


2.365 


2.226 


Wc 


9.7(1) 


7.98(02) 


5.77(02) 


V 


2.27(06) 


2.29(08) 


2.82(05) 


21 


2.04(04) 


1.59(01) 


0.91(01) 


<9> 


2.00 


2.82 


7.99 


var g 


0.84 


1.24 


2.53 




0.58 


0.96 


2.04 


var log g 


0.31 


0.22 


0.56 



Table 1: Critical parameters of bifractals A, B and C 
(Figure 1). d s is spectral dimension |2(J. Critical dis- 
order and critical exponent were found from fit of the nu- 
merical data by a simple one-parameter scaling formula. 
Accuracy of critical parameters was estimated by com- 
parison of numerical fits for system size L = 4 32 and 
L = 8-4-32 and for various intervals of W. Limitting values 
of the moments of conductance were estimated by the fit 
X L = + const/L with L = 8, 16, 32 and 64 (figure 2). 



lattice 


3s 


3h 


3t 


4 


d a 


3 


3 


3 


4 


Wc 


16.5 


13.5(2) 


19.9(2) 


34.3(2) 


V 


1.54(03) 


1.58(04) 


1.53(03) 


1.1(1) 


Zl 


3.45(01) 


2.7(3) 


4.2(2) 


5.4(1) 


(9) 


0.57 


0.96 


0.30 


0.27 


var g 


0.17 


0.58 


0.19 


0.17 


(log a) 


-1.93 


-0.44 


-2.34 


-2.39 


var log g 


1.76 


1.00 


3.02 


2.61 



Table 2: Critical parameters of three dimensional systems 
with square (3s) |^|, H|| , honeycomb (3h) and triangular 
(3t) lattice in xy plane and of 4D hyper-cube. Periodic 
boundary conditions (b.c.) were used for Q1D simulations. 
Fixed b.c. were used for studies of conductance statistics. 



